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Abstract 

The continuation of resonant periodic orbits from the restricted to the general 
three body problem is studied in a systematic way. Starting from the Keplerian un- 
perturbed system we obtain the resonant families of the circular restricted problem. 
Then we find all the families of the resonant elliptic restricted three body problem 
which bifurcate from the circular model. All these families are continued to the 
general three body problem, and in this way we can obtain a global picture of all 
the families of periodic orbits of a two-planet resonant system. We consider planar 
motion only. We show that the continuation follows a scheme proposed by Bozis 
and Hadjidemetriou (1976) for symmetric orbits. Our study includes also asym- 
metric periodic orbits, which exist in cases of external resonances. The families 
formed by passing from the restricted to the general problem are continued within 
the framework of the general problem by varying the planetary mass ratio p. We 
obtain bifurcations which are caused either due to collisions of the families in the 
space of initial conditions or due to the vanishing of bifurcation points. Our study 
refers to the whole range of planetary mass ratio values {p € (0, oo)) and, therefore 
we include the passage from external to internal resonances. In the present work, 
our numerical study includes the case of the 2/1 and 1/2 resonance. The same 
method can be used to study all other planetary resonances. 

Keywords Three body problem (TBP), Resonances, Continuation of periodic orbits, 
Bifurcations. 



1 Introduction 

A good model to study the motion of three celestial bodies considered as point masses, e.g. 
a triple stellar system or a planetary system with two planets, is the famous three-body 
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problem (TBP), whose study goes back to Poincare. In the present study we consider 
the case where only one of the three bodies is the more massive one, and the other two 
bodies have much smaller masses. This model is useful in the study of the motion of small 
bodies in a planetary system (for example asteroids and comets in our Solar System) or 
planetary systems with two planets. 

The simplest model is the circular restricted TBP. Although much work has been 
carried out for this model, (see e.g. Bruno, 1994; Henon, 1997), new interesting re- 
sults continue to appear in the literature (Maciejewski and Rybicki 2004; Papadakis and 
Goudas, 2006; Bruno and Varin 2006,2007). In this model we consider two bodies with 
non zero mass, called primaries, for example the Sun and Jupiter, moving in circular 
orbits around their common center of mass, and a third body with negligible mass, for 
example an asteroid, moves under the gravitational attraction of the two primaries. A 
more realistic model is the elliptic restricted TBP, where the two primaries move in el- 
liptic orbits, with a finite eccentricity. However, an important aspect of the dynamics is 
missing in these two models. It is the gravitational interaction between the small body 
and the two primaries, which is not taken into account in the restricted models. When 
we introduce this gravitational interaction, we have a more realistic model, the general 
TBP. Within this framework we study a system consisting of a Sun and two small bodies, 
which we call planets. 

In the study of a dynamical system, the topology of its phase space plays a crucial role. 
The topology is determined by the position and the stability properties of the periodic 
orbits, or equivalently, of the fixed points of the Poincare map on a surface of section. 
This makes clear the importance of the knowledge of the families of periodic orbits in a 
dynamical system. Particularly, in a planetary system, many families of periodic orbits are 
associated with resonances, which are mean motion resonances between the two planets. 
Since in our study of planetary systems only one body is the more massive one (the 
sun), a good method is to start from the simplest model, which is the circular restricted 
problem and find all the basic families of periodic orbits. Then we extend the model 
to the elliptic restricted model, and find all the families of resonant periodic orbits that 
bifurcate from the circular to the elliptic model. Finally, we give mass to the massless 
body and continue all these families to the model of the general problem. This is the 
method that we shall use in the present study to describe the topology of the phase 
space. The existence of periodic orbits of the planetary type in the general TBP, as 
a continuation from the restricted problem, has been studied by Hadjidemetriou (1975, 
1976) and recently this method found a fruitful field of applicability in the dynamics 
of resonant extrasolar systems (e.g. Rivera and Lissauer,2001; Hadjidemetriou, 2002; 
Ji et al, 2003; Haghighipour et al, 2003; Ferraz-Mello et al., 2003, 2005; Psychoyos 
and Hadjidemetriou, 2005; Voyatzis and Hadjidemetriou, 2005, 2006; Voyatzis, 2008). 
An alternative method for determining resonant periodic orbits is the computation of 
stationary solutions of an averaged model (Beauge et al. 2003; Michtchenko at al. 2006). 

It is interesting to mention at this point that the situation is not the same in all 
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resonances. It may happen that in some resonances there are not new families in the 
elliptic model, bifurcating from the circular model, or may exist several resonant families 
bifurcating from the circular model. This has, evidently, important consequences on 
the topology of the phase space and differentiates the behavior of the model in different 
resonances (see Tsiganis et al. 2002a, 2002b). Although in this study we present numerical 
computations of the 2/1 (or 1/2) resonance, our approach and results have a more general 
applicability. 

In the next section we discuss briefly known aspects on the families of periodic orbits in 
the circular unperturbed and restricted problem. However, these issues are fundamental 
for continuing our study in the elliptic restricted and in the general problem. In section 
[3] we present the continuation of resonant periodic orbits from the circular to the elliptic 
model and then, in section HJ we consider the continuation in the general problem. In 
section we study the bifurcation of families of periodic orbits within the framework of 
the general problem and finally, in section |6j we discuss the generality of our results and 
conclude. 

2 The circular restricted problem 

Consider a body S (Sun) with mass m and a second body J (Jupiter) with mass mi, 
which describe circular orbits around their common center of mass. We define a rotating 
frame of reference xOy, whose x-axis is the line SJ, the origin is at their center of mass 
and the xy plane is the orbital plane of the circular motion of the these two bodies. The 
circular restricted problem describes the motion of a massless body A in the rotating 
frame, which moves under the gravitational attraction of S and J. In our computations 
we consider the normalization of units m + mi = 1, G = 1 and n' = 1, where G is the 
gravity constant and n' the mean motion of J. 

2.1 The unperturbed case (mi = 0) 

If J is assumed to be a massless body too, then in the inertial frame the body A can move 
in a Keplerian ellipsis of semimajor axis a, eccentricity e and mean motion n ~ a~ 3//2 . 
Circular orbits (e = 0) with radious r = r = a are also included. The above mentioned 
orbits are periodic in the inertial frame. 

In the rotating frame, the Hamiltonian that describes the unperturbed motion of A, 
in polar coordinates, r, 0, is 
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The momenta are p r = r and p^ = r 2 (<f) + n')=constant. In the rotating frame periodic 
orbits also exist and are of two kinds: 
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Figure 1: (a) The unperturbed case where Jupiter mass is zero (mi = 0). The family of 
circular periodic orbits and the the two resonant families at the 2/1 and 1/2 resonance that 
bifurcate from the circular family, (b) The bifurcation to the two 2/1 resonant elliptic 
families (inner resonance), when the mass of Jupiter is non zero. The same topology 
appears at the bifurcation point at the 1/2 resonance (outer resonance). 



Circular orbits. In the rotating frame there exist circular orbits of the body A with an 
arbitrary radius r = r and are evidently symmetric with respect to the rotating x-axis. 
Consequently, a family of circular periodic orbits exists along which the radius ro = a or 
the frequency n = (GmoY^a^ 3 ^ 2 varies. This family is represented by a smooth curve 
in the space Xq — h, where Xq is the initial condition of the orbit, which is Xq = r or, 
equivalently, Xq = a, and h denotes the value of the energy integral H . The family is 
given by the equation 

- - - nV Gm x = h, (2) 
2x 

Elliptic orbits. An elliptic orbit of the small body A in the inertial frame is periodic in 
the rotating frame only if it is resonant, i.e. n/n' = p/q, with p,q integers. If ro = a\ 
is the semimajor axis (radious) of Jupiter then a n/n' = p/q resonance corresponds for 
the massless body A to a semimajor axis a = a p / q = (n/n')~ 2 ^ 3 and is periodic for any 
eccentricity e. Thus & family of resonant elliptic periodic orbits exists, with the eccentricity 
as a parameter along the family. The ratio n/n' is almost constant along the family. There 
is however another parameter, defining the orientation of the elliptic orbit, which is the 
angle u of the line of apsides with a fixed direction. If this fixed direction is defined by a 
conjunction of the three bodies then an elliptic periodic orbit is symmetric when u = 
or u = 7i. 

In the space x — h, an elliptic family is represented by a smooth curve, given by 
the equation (J2]) where now xq = a p / q (l — e) is the pericenter distance. Note that this 
presentation is not unique: a point on the elliptic family represents all the elliptic resonant 
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orbits with the same eccentricity, but arbitrary orientation u. An elliptic periodic orbit 
in the rotating frame is also periodic in the inertial frame. Along the circular family the 
value of the semimajor axis varies and, consequently, the ratio n/n' varies and passes 
through resonant values a = a p / q . It is at these points that we have a bifurcation to an 
elliptic family. 

The circular family and the 2/1 and 1/2 resonant families of periodic orbits of the 
unperturbed problem are presented in Fig{Th. In the normalization we are using, the semi 
major axes of the 2/1 family are all a < 1, called internal resonance and the semimajor 
axes of the 1/2 resonance are all a > 1, called external resonance. The bifurcation points 
are r = a 2 /i = 0.6297 and r = ai/ 2 = 1.5874, respectively. At the bifurcation points, 
the tangent to the above resonant elliptic families is parallel to the x-axis. 

Each family of the internal and the external resonance, is divided into two parts by 
the corresponding bifurcation point. One part corresponds to position of the small body 
at perihelion and the other at aphelion. In particular, for the internal family, the part 
x < 0.6297 corresponds to perihelion (family Hi) and the part x > 0.6297 corresponds 
to aphelion (family !{). For the external family, the part x < 1.5874 corresponds to 
perihelion (family I e ) and the part x > 1.5874 corresponds to aphelion (family II e ). Note 
that on the aphelion part of the inner family and on the perihelion part of the outer 
family, there is a collision orbit with Jupiter. Evidently, all the circular and the elliptic 
orbits are stable, as they are Keplerian orbits. 

2.2 The perturbed case: Non zero mass of Jupiter 

Let us now assume that the mass of Jupiter is non zero i.e. m\ = fi ^ and mo — 1 — fi. 
The resonant elliptic families at the 2/1 and 1/2 resonance are continued to fi > 0, but 
a gap appears at the bifurcation point of the unperturbed families . The topology at the 
bifurcation point of the resonant families and the gap that appears, is shown in Figure [2b. 
There are two families of elliptic orbits in each resonance, one corresponding to the case 
where the small body is at perihelion and the other to the case where the small body is at 
aphelion. All these orbits are symmetric periodic orbits, i.e. out of the infinite set of all 
symmetric and asymmetric orbits of the unperturbed problem, only two orbits survive for 
/i > 0, both symmetric. One of them is stable and the other unstable, as a consequence of 
the Poincare-Birkhoff fixed point theorem, but along the family the stability may change. 

In Figure [2^ we present the families of periodic orbits for the 2/1 internal resonance, 
for fi = 0.001. The gap at the bifurcation point that we showed in Figure |2b is presented 
in this Figure by a small gap at the point e\ ~ ix rs 0.6297). The families are presented 
in a plane where the horizontal axis indicates the x-coordinate of the massless body in 
the rotating frame and the vertical axis corresponds to the osculating eccentricity of the 
nearly Keplerian orbit of the massless body. Correspondingly to the unperturbed case, 
there are two families, family 1^, where the small body is at aphelion and family 
where the small body is at perihelion. Family La consists of two parts: L ia ) and 
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Figure 2: Resonant families of elliptic periodic orbits of the circular restricted problem, 
a) the case of internal 2/1 resonance, b) the case of external 1/2 resonance. Thin or bold 
curves indicate unstable or stable orbits, respectively. The symbol " x" indicates a region 
of close encounters between Jupiter and the small body. Note that we denote with e\ the 
eccentricity of the small body in the inner orbits and with €2 the eccentricity of the small 
body in the outer orbits. 

separated by a collision (see Figure |2h). The first part consists of unstable orbits (thin 
curve) while the second part consists of stable orbits (thick curve). The orbits of the 
family IIu) are all stable. 

In the case of the external 1/2 resonance, the corresponding families are shown in 
Figure [2b. These families arise from the 1/2 elliptic family of Figure |2k when the last 
one is continued to fi >0. As in the internal resonance 2/1, a small gap appears at the 
point C2 = 0, x = 1.5874. There are two families, family Ir e \, where the small body is 
at perihelion and family II( e ), where the small body is at aphelion. As in the internal 
resonance 2/1, the family L e \ consists of an unstable part (/( ea )) an d a stable part (ir e &)), 
which are separated by a collision (see Figure Ek>). The family IIr e \ starts and ends 
with stable orbits, but along the family the stability changes and an unstable part exists 
between the critical points B x c and Bq. From each one of these critical points there 
bifurcates a family of asymmetric periodic orbits. It turns out that these two asymmetric 
families coincide to a single asymmetric family A, which starts from the point B\, and 
ends to the point B^,. We have found that this family is whole stable for \x < 5.210~ 3 . 
We remark that the asymmetric family exists only for the external resonances of the form 
1/q, called asymmetric resonances (Beauge, 1994; Voyatzis et al., 2005). 
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3 The elliptic restricted model 



Along the resonant elliptic families of the circular restricted problem the period varies. In 
the unperturbed case (/x = 0) the period is exactly equal to 2n along the family, for the 
normalization we are using. When \i > 0, the period along the elliptic resonant families 
varies, but is close to the value 2n. If it happens that for a particular orbit of the family 
the period is exactly equal to T = 2tc, this point is a bifurcation point for a family of 
periodic orbits of the elliptic restricted TBP, along which the eccentricity of the second 
primary (Jupiter) varies. The same is true in the general case where the period of a 
periodic orbit on a family of the circular restricted problem is a multiple of tt (Broucke, 
1969a,b). 

Concerning the stability, we remark that the monodromy matrix of the variational 
equations has two pairs of eigenvalues (Ai,A2) and (A3, A4). In the circular model one 
pair, (Ai, A2) is the unit pair, Ai = A2 = 1, because of the existence of the energy integral. 
The other pair, (A3,A4), may lie on the unit circle in the complex plane (A3 : 4 = e ±J<7i ), 
corresponding to stability, or may be on the real axis (A3 = I/A4 G R) corresponding 
to instability. In the elliptic model there is no energy integral, so we have four different 
cases: (i) stable orbits when all eigenvalues are on the unit circle (ii) simply unstable 
orbits when one pair of eigenvalues is on the unit circle and one on the real axis (iii) 
doubly unstable orbits when both pairs of eigenvalues are on the real axis and (iv) complex 
instability, where all eigenvalues are outside the unit circle, not on the real axis, arranged 
in reciprocal pairs and complex conjugate pairs Ai 2 = Re ±l ^, A 3i4 = R~ x e ±%< ^ (see Broucke 
1969a,b). We also can define the stability indices b\ = Ai + A2 and 62 = A3 + A4 (see 
Hadjidemetriou, 2006). A periodic orbit is stable if 

< 2, Vi = 1,2. 

If one of the indices bi does not satisfy the stability condition, the periodic orbit is simply 
unstable. If both indices do not satisfy the stability condition the periodic orbit is doubly 
unstable. Complex instability is not present in the particular model. 

In the case of 2/1 and 1/2 resonant families, the critical points, i.e. periodic orbits 
with period exactly equal to 27r, are indicated on the families IIu\ and H( e ) in Figure 
[2] by the points B l T , i = 0, 1, 2, ... 6. Note that for the internal resonance 2/1, there is 
only one bifurcation point, Bj,, on the stable family Hu). For the external resonance 1/2 
there are five bifurcation points. Three of them, B)p, B\ and B\ belong to the symmetric 
family 11(e) an d two more points, B\, B\ belong to the asymmetric family A. There are 
no bifurcation points on the families Iu\ while on the family there is one bifurcation 
point, B\ at high eccentricity value. 

Let us consider a critical periodic orbit B l T of the circular problem. The continuation 
of a family of periodic orbits of the elliptic model, which bifurcates from a point B % Ty 
is obtained by increasing the eccentricity of Jupiter, starting from the zero value and 
keeping its semimajor axis equal to unity, 0! = 1, so that n' — 1. This is a family 
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Figure 3: Families of 2/1 resonant periodic orbits in the elliptic restricted problem. The 
critical orbits B l Tl i = 0,1,2,4,5 of the circular problem are indicated. Bold, thin and 
dotted face curves indicate stable, unstable and doubly unstable orbits respectively. 

of resonant periodic orbits, along which the eccentricity of Jupiter increases. There are 
two possibilities (Hadjidemetriou, 1993): the elliptic orbit of Jupiter corresponds (i) to 
perihelion at t = or (ii) to aphelion. So, there are two families of resonant periodic 
orbits of the elliptic model which bifurcate from each point B l T . In one family, denoted by 
E p , Jupiter is at perihelion and in the other family, denoted by E a , Jupiter is at aphelion. 
Additionally, as we will show in the following, families of asymmetric periodic orbits also 
exist for the elliptic problem and can be generated in three different ways. We will denote 
the asymmetric orbits by E A . 

The resonant families of the elliptic model will be presented in the space of the eccen- 
tricities of Jupiter and the small body. There are two bodies involved (Jupiter and the 
small body) and in all the following we will call e\ the eccentricity of the inner of these 
bodies and e2 the eccentricity of the outer body, irrespectively of whether the correspond- 
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ing body is Jupiter or the massless body. We use this notation in order to have a direct 
comparison with the families of the general TBP presented in the following sections. 

Let us start from the family of the circular model, and present the family of the 
elliptic model that bifurcates from the orbit Bj, (Figure Ek.), in the space e\-e 2 . Now e\ 
denotes the eccentricity of the small body (inner body) and e 2 denotes the eccentricity of 
Jupiter (outer body). The Family 11^ of the 2/1 resonant family of the circular model is 
located in Figure E] on the axis e 2 = (circular orbit of Jupiter). On this axis we present 
the point B^ and show the two families of the elliptic model that bifurcate from this point. 
As we mentioned before, there are two families, Eo a and Eq p , corresponding to perihelion 
and aphelion of Jupiter, respectively. Both families start from the eccentricity e\ = 0.75 
of the small body and the eccentricity e 2 of Jupiter starts from zero and increases along 
the family. 

We come now to the families of the elliptic model that bifurcate from the 1/2 external 
resonant family II( e ) of the circular model (Figure [2b). In this case, e\ is the eccentricity of 
Jupiter (inner body) and e 2 is the eccentricity of the small body (outer body). The family 
J/( e ) is located in the Figure [3] on the axis t\ = 0. Let us consider first the bifurcation 
from the symmetric families of the circular problem. We study the two critical points B\ 
and B\ (the point B\ corresponds to very high eccentricities and we do not study it here). 
From the point B\ there bifurcate two resonant families of symmetric periodic orbits of 
the elliptic model, the family E\ p , which starts having stable orbits, and the family E\ a 
where the orbits are doubly unstable. From the critical point B\ there also bifurcate two 
families of resonant symmetric periodic orbits of the elliptic model, the families E 2a and 
E 2p ; both are unstable. It turns out that the family Ei a that bifurcates from the point B\ 
and the family E 2a that bifurcates from the point B\ coincide and form a single family 
that starts from B\ and ends to B\. 

We come next to the asymmetric families of the elliptic model. There exist two critical 
points on the asymmetric family A in FigfJJ the points B\, and B\. As in the previous 
cases, from each one of these points we have a bifurcation of two resonant families of the 
elliptic model, which are asymmetric. From the point B\ we have the bifurcation of the 
family Ef x (stable) and the family Ef 2 (unstable). This latter family has a complicated 
form, its stability changes three times and ends at the point B\. So it seems that from this 
latter point there bifurcate, in addition to the two symmetric families, one more asym- 
metric family. From the asymmetric point B\ there bifurcate two asymmetric families, 
and E^ 2 . The family E^ x is stable but family E^ 2 cannot be numerically continued 
due its strong instability and, thus, is not presented in Figj2j 

Along the resonant families of the elliptic model it may also happen to exist bifur- 
cation points. Such a case is with the family E ip , where the critical point B l E appears. 
This happens because the stability type changes at this point. From this point we have 
a bifurcation of the family A^ p of asymmetric periodic orbits. Voyatzis and Kotoulas 
(2005) showed that many families of the external resonances have such critical points and 
conjectured the bifurcation of asymmetric orbits. 
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4 From the restricted to the general problem 



4.1 The rotating frame and periodicity conditions 

Let us consider three bodies with non zero masses, with one of them much more massive 
than the other two. The more massive body with mass mo is the Sun, S, and the two small 
bodies, Pi with mass mi and P 2 with mass m 2 , will be called planets. In the following, 
P 1 will be the inner planet and P 2 the outer planet. If mi 7^ and m 2 = we have 
the restricted model where Pi is the corresponding Jupiter and P 2 the small body, which 
moves initially in an outer orbit. If the Keplerian orbit of the bodies S and Pi is circular, 
we have the circular restricted problem and if it is elliptic, we have the elliptic restricted 
model. If mi = and m 2 7^ then P 2 plays the role of Jupiter and the massless body Pi 
evolves initially in an inner orbit. 

In order to study the continuation of the periodic orbits from the restricted to the 
general problem we define a rotating frame of reference xOy, whose x-axis is the line 
S — Pi, with the center of mass of these two bodies at the origin O and the y-axis is in 
the orbital plane of the three bodies. In this rotating frame the body Pi is always on the 
x-axis and P 2 moves in the xOy plane. We have four degrees of freedom and we use as 
coordinates the position x\ of Pi, the coordinates x 2 , y 2 of P 2 and the angle 9 between 
the x axis and a fixed direction in the inertial frame. 

The Lagrangian in the above coordinates is (Hadjidemetriou 1975) 



mi 



+V, 



_ ^.2 + r 2^ 2) + rrn U 2 + 2 + §2 ^ 2 + y2) + ^ _ A 
mo ' m V / 

(3) 

where m = m + mi + m 2 is the total mass of the system, /i — m\j (m + mi), r — x±/ fj, 
is the distance between the Sun and Pi and 

Gm mi Gm\m 2 Gm m 2 



r v /((i_ /i ) r _ a;2 )2 + y 2 y/(fj,r + x 2 ) 2 + yZ' 

Note, that 9 is ignorable and consequently the angular momentum L = dC/d9 is constant. 
We can use this angular momentum to reduce the number of degrees of freedom from four 
to three. In the reduced Lagrangian, called the Ruthian, the angular momentum constant 
L is a fixed parameter. For this reason we take as normalizing conditions to fix the units 
of mass, length and time the conditions 

m = 1, G — 1, L — constant. 

Periodic orbits of period T exist in the rotating frame and their initial conditions must 
satisfy the periodic conditions (Voyatzis and Hadjidemetriou, 2005): 

Xl {T) = Xl {0) = 0, x 1 (T) = x 1 (0), 
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x 2 (T) =x 2 (0), y 2 (T)=y 2 (0), 
x 2 (T)=x 2 (0), y 2 {T)=y 2 {0). 



Thus, a periodic orbit is represented by a point in the 5- dimensional space 

n = {(x 1 (0),x 2 (0),t/ 2 (0),x 2 (0), y 2 (0))}. (4) 

When a periodic orbits is symmetric, i.e. it is invariant under the fundamental symmetry 
£ : (t,x,y) — > (—t,x,—y) (Henon, 1997; Voyatzis and Hadjidemetriou, 2005), we can 
always take as initial conditions y 2 (0) = and x 2 (0) = and the dimension of the space 
of initial conditions II is reduced to three. In the planetary problem, the orbits of Pi and 
P 2 are almost Keplerian and we can present the families of periodic orbits in the projection 
plane of planetary eccentricities e\ — e 2 , which correspond to the initial conditions. 

The monodromy matrix of the variational equations has now three pairs of eigenvalues 
(Ai,A 2 ), (A3, A 4 ) and (A 5 , A 6 ). Due to the existence of the energy integral it is A 5 = A 6 = 1 
and the stability of the periodic orbits is defined by the first two pairs as in the elliptic 
restricted problem (see section [3]). 

4.2 The continuation from the restricted to the general problem 

In general, the periodic orbits of the restricted problem are continued to the general 
problem. Particularly, it is proved by Hadjidemetriou (1975) that a periodic orbit of 
period T of the circular restricted problem is continued to the general problem with 
the same period, by increasing the mass (e.g. m 2 ) of the initially massless body. The 
continuation is not possible only in the case where the period is a multiple of 27r (the period 
of the primaries). Such continuation forms monoparametric families of periodic orbits with 
parameter the mass of the small body (m 2 ), provided that the masses of the other two 
bodies are fixed. If we keep all masses fixed (and non zero), we obtain a monoparametric 
family of periodic orbits of the general problem, of the planetary type, along which the 
elements of the two planetary orbits vary. This family is represented by a smooth curve 
in the space IT of initial conditions. The periodic orbits of the elliptic restricted problem 
are also continued to the general problem (Hadjidemetriou and Christides, 1975). 

The evolution of the characteristic curves of the families, as m 2 increases, is studied 
by Bozis and Hadjidemetriou (1976). In the present paper we study the continuation 
using a slightly different approach. Provided that the planetary masses are small with 
respect to the mass of Sun, i.e. mi <C m and m 2 <C m , it is found that the families of 
resonant periodic orbits in the space n depends on the ratio p = m 2 /m 1 of the planetary 
masses rather than by their actual values (Beauge et al., 2003). Thus, we can use p as 
a continuation parameter for a monoparametric family of periodic orbits. If we add one 
more dimension to the space of initial conditions n in order to assign the value of p, we 
obtain an extended space n'. In this extended space we can form characteristic surfaces 
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Figure 4: The schematic generation of families of periodic orbits in the general problem 
from the restricted problem. Dashed curves indicate the families of the restricted problem 
and the solid curves indicates the families continued in the general problem (see text for 
detailed description). 

of two-parametric families. In the following we will present families of periodic orbits 
considering sections of II' defined by fixing p to a constant value. 

In Fig. H] it is shown schematically the continuation of the families F and E of the 
circular and the elliptic restricted problem, respectively, with e.g. p = 0. The point 
A of the family F of periodic orbits of the circular model corresponds to an orbit with 
period T = 2k-K, where k is an integer. At this point the families E p and E a bifurcate 
to the elliptic restricted problem. By giving mass to the small body, any orbit of F with 
T ^ 2kir, e.g. the orbit B, is continued to the general problem with parameter the mass 
ratio p and the monoparametric family g F is formed. A family g F exist for any initial 
orbit of F, except for the orbit A. If we start now from a point on the family g F and 
keep p fixed, we obtain a monoparametric family c. The same continuation holds for the 
periodic orbits of the families E p and E a of the elliptic restricted problem and for a fixed 
p 7^ we get the characteristic curve d. The individual parts of c and d join smoothly 
to each other forming the families G p for the general problem and for any fixed mass 
ratio p C 1. The periodic orbit A is a singularity for the continuation resulting in the 
formation of a gap between the left and right characteristic curves cd. 

4.3 Families of 2/1 and 1/2 orbits continued from the restricted 
problem to the general problem 

Concerning the continuation of the resonant families, which are presented in sections 12.21 
and [31 to the general problem, singularities are expected at the bifurcation points (see 
Fig. [3} for both internal (2/1) and external (1/2) resonances, where is T = 2kn. In the 
following, we study four different cases of continuation presented in Fig. [5j We remind 
that the families of the circular restricted problem are continued as curves close to the 
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vertical and the horizontal axes in the plane e\ — e^, for the external (ei = 0) and the 
internal resonances [e^ = 0), respectively. In our computations we vary p by changing the 
value of one of the planetary masses and such that max(mi,m2) = 1CT 3 . 

In the panel (a) of Fig. [5] we show the continuation close to the bifurcation point Bj>, 
which belongs to the symmetric family II ( e ) of the circular problem and is the starting 
point for the families E\ p and Ei a of the elliptic problem. By giving mass to the small 
body (i.e. p ^ 0) the part (OB^) of the family II( e ) and the family E\ p are continued 
and they join smoothly forming the symmetric family Si. The change of stability that 
is shown in family H( e ) at the point B\, is also obtained in the family Si. Namely, 
the bifurcation point B\, (periodic orbit of critical stability) continues as a bifurcation 
point B' c in the general problem too. The asymmetric family A of the circular restricted 
problem, which bifurcates from B^, is continued smoothly to the general problem as an 
asymmetric family which bifurcates from the B' c . Similarly to the formation of Si, the 
unstable part of the family H( e ) above the point B\ and the family Ei a , which is doubly 
unstable, are continued and form the family S2. Note that the different stability types of 
the two families results to a periodic orbit of critical stability in the family S2. 

In family E\ p there is a change in stability at a point Be, where the asymmetric family 
E^ bifurcates for the elliptic problem. As p takes a positive value, B E continues as a 
critical point in the general problem and the asymmetric family A3 bifurcates from it 
(Fig. [5b). Actually, A3 can be considered as the continuation of the family Ef. No 
singularities are obtained in this case. Such a type of continuation explains the existence 
of the asymmetric planetary corotations found by Michtchenko et al. (2006) at the 3/2 
resonance, which in the elliptic restricted problem shows an asymmetric family similar to 
Ei p . We note that asymmetric periodic orbits were known only for resonances p/q with 
q = l. 

As it is mentioned above, the asymmetric family A continues smoothly to the general 
problem (as family A\) starting from its bifurcation point B\, (see Fig. . However, 
the family A contains the critical point B^ (see Fig. [2]), which admits a singularity for 
the continuation to the general problem. In the panel (c) of Fig. [5j the gap formed in 
the above singularity is shown. The asymmetric family E^_, which bifurcates from B^, 
continues and completes the asymmetric family A\ of the general problem. Additionally, 
the part of family A located between the bifurcation points B\ and B^, the family Ef 2 and 
the family E^ continue in the general problem, join smoothly and form the asymmetric 
family A 2 . Again we obtain a change of stability along A 2 due to the different types of 
stability on the families of the restricted problem (family A is stable while family E^ 2 is 
unstable). 

In the case of the internal resonance, the continuation to the general problem shows 
the same characteristics as in the external resonances. We remark that in this case only 
symmetric periodic orbits appear. In Fig. [5H we show the formation of the symmetric 
family S[ of the general problem, after the continuation of the families 11^ and Eq p of 
the restricted problem. We note that since p = 00 in this case, the continuation maybe 
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assumed as varying 1/ p. 



5 Continuation within the framework of the general 
problem 

In the previous section we considered the continuation of families of periodic orbits when 
we pass from the restricted to the general problem. The situation, which was described, 
is referred to small values of the mass of one of the planets, namely p <C 1 or p ^> 1. Now 
we consider the case of the external resonance and its families generated for p > 0. As 
the parameter p increases, the families evolve and new bifurcations and structure changes 
are possible. We follow the evolution of the 2/1 resonant families, which are presented 
in the previous section, and we show in the following how the new structures are formed 
by increasing p. We remind that as p — > oo we approach the internal resonance of the 
restricted problem. 

Figure shows the asymmetric families A 2 and A 3 for p = 0.2. Up to this value not 
any structural change of the families occur. In both families the stability changes and, 
consequently, we distinguish in each family two parts, parts a and b for A 2 and c and d for 
A%. For the critical value p = p\ ~ 0.275 the two families collide at a point in the space of 
initial conditions II (Fig. [6b). This point corresponds to the orbit of critical stability and 
is the intersection point of the parts a-d. For p ^> pi (see Fig. [6b) the part a of family A 2 
and the part d of family A 3 join together and form the family A23. Similarly, the parts 
c and b of the families A 3 and A 2 , respectively, form the family A 32 . Note that the new 
families are separate and their intersection in the plane e\ — e 2 is due to the projection. 

The evolution of the family A\, as p increases, take place smoothly without structural 
changes up to p = p 2 ~ 0.37. In Figure^, which corresponds to p = 0.3, it is shown that 
the family Ai has come close to the family A 32 , which is generated after the bifurcation 
at p = p\. At p = p 2 the two families collide and, similarly to the previous case, for 
p > p 2 we obtain two new families, the family A4 and the family A123 (Fig. [7b). In this 
case only the family A% 2 has an orbit of critical stability, which separates the family in 
a stable (a) and in an unstable part (b). The family A\ is whole stable and there is no 
clear border between its parts c and d. After the bifurcation, an orbit of critical stability 
is shown only in the new family A^, while the family Ai 23 is whole stable and starts and 
ends at bifurcation points of the symmetric family Si. Note that these bifurcation points 
originate to the bifurcation points B\, and Be of the circular and the elliptic, respectively, 
restricted problem. 

Now we restrict our study to the evolution of the family Ai 23 for p > p 2 . As it is shown 
in Fig. [Sk, as p increases, the ending points of the family move on along the family Si in 
opposite direction and the family shrinks and, finally disappears at p = pz w 1.034. In 
Fig. [8b we present the above transition by considering the stability indices 61,62 for the 
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orbits along the family S\ (see Hadjidemetriou, 2005). The horizontal axis indicates the 
eccentricity of the periodic orbits along the family Si and the value of the corresponding 
stability index b 2 is presented on the vertical axis. For all orbits it holds |6i| < 2 and, 
thus, the condition \b^\ < 2 is sufficient and necessary for linear stability. In Figure [Hb we 
obtain that unstable orbits exist for p < p 3 in a part of the family S± defined by the e\ 
interval where b 2 < — 2. As p increases, the curve of 62 values is raised continually and for 
p > p 3 is located above the value b 2 = —2. Therefore, the unstable part of Si disappears 
and, consequently, the family Ai 23 disappears too. 

Actually for p > p 3 , the planet Pi becomes the small body and we pass to the case 
of the internal resonance. The remaining family Si evolves smoothly, as p increases, and 
should be assumed as a family S[ of the internal resonance, which approaches the families 
1 1 (a and E 0p of the restricted problem as p — > 00 (see sections I2T2151 and Fig. |SH). 

The above described bifurcation scheme explains completely the origin and the struc- 
ture of corotations found by Beauge et al. (2006). Following the other asymmetric fami- 
lies, which are coming from the restricted problem, we obtain new "collision-bifurcations" 
as p increases and new structures of characteristic curves are formed. We have found e.g. 
that the families and A 23 collide for p 0.45 and generate the miscellaneous family 
found by Voyatzis and Hadjidemetriou (2005) and called ll A 2 " there in. 

6 Discussion and Conclusions 

The method of continuation of periodic orbits has been applied in the present work in 
order to study the generation and the structure of families of periodic orbits in the general 
TBP, starting from the restricted problem. We considered the planar TBP of planetary 
type, referred to a rotating frame and studied resonant planetary motion. In particular, 
we studied the 2/1 resonant periodic orbits, but the method we used can be applied to 
all other resonances. We considered planar motion only. 

We started from the circular restricted TBP and computed the basic families of 2/1 
resonant periodic orbits, both for the inner orbits (inside Jupiter) and the outer orbits 
(outside Jupiter). Along these families the resonance is almost constant, but the eccen- 
tricity of the small body increases and may take high values. The basic families are 
symmetric with respect to the rotating x-axis, but asymmetric families also exist in cases 
of resonances of the form \jq. On these families, we found the periodic orbits with period 
equal to 2tt which, in the normalization we are using, are the bifurcation points to resonant 
2/1 (or 1/2) families of periodic orbits of the elliptic restricted model, along which the 
eccentricity of Jupiter varies, starting from zero values. Two such families bifurcate from 
each of these critical points. In particular, we have the following cases: a) an asymmetric 
periodic orbit of period 2ir of the circular restricted problem is continued to the elliptic 
restricted problem as asymmetric also, b) a symmetric periodic orbit of period 2it of the 
circular restricted problem can be continued to the elliptic problem as an asymmetric one 
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- this exceptional case is verified only up to the accuracy of the numerical computations 
and c) a symmetric periodic orbit of the elliptic restricted problem, which is of critical 
stability, can be continued to the elliptic problem as an asymmetric one. The last case 
indicates the existence of asymmetric periodic orbits in the elliptic problem and in par- 
ticular to resonances which are not necessarily of the form 1/q. In this way we obtain a 
clear picture of all the resonant families of the restricted problem, and in this way we can 
make a complete study of the resonant families of the general problem, by giving mass to 
the massless body of the restricted model. 

The continuation of the families of periodic orbits of the restricted problem to the 
general problem follows the scenario described in the paper of Bozis and Hadjidemetriou 
(1976). Particularly, the families of the general problem originate from two families of the 
restricted problem, one of the circular problem and one of the elliptic problem. There is 
no essential difference in the continuation between symmetric and asymmetric periodic 
orbits. The asymmetric orbits, which bifurcate from symmetric orbits of the elliptic 
problem, are continued smoothly in the general problem. The stability of periodic orbits 
is preserved after the continuation from the restricted to the general problem. In the paper 
mentioned above, a case of continuation without the formation of a gap is described. We 
found this case only at bifurcation points of the families of the elliptic restricted problem 
(see FigfSb). 

After the continuation of the periodic orbits of the restricted problem to the general 
problem, by giving a very small mass to the massless body, we studied how these families 
evolve when we increase the mass of the small body. Particularly we studied such an 
evolution by considering as a parameter the planetary mass ratio p = 7712/1711, keeping 
mi -C l,i = 1,2. In our numerical computations we considered max(mi, 777.2 ) = 10" 3 and 
we varied p by varying only the mass of the smaller body. Starting from p = (external 
resonances of the restricted problem) and increasing its value, we found that the charac- 
teristic curves of two different families can collide in the space of initial conditions. At 
these points a bifurcation takes place (collision-bifurcation) causing a topological change 
in the structure of the colliding families and the formation of new families. As p — > 00 
we approximate the families of the internal resonances of the restricted problem. Since 
asymmetric periodic orbits are not known for the internal resonances, it would be inter- 
esting for a future work a study of the evolution and the bifurcation of all asymmetric 
families up to values of p > 1 and, furthermore, as p — > 00. 

In general, the continuation and the evolution the families of the external resonances 
(p < 1) in the general TBP show more rich dynamics (number and structures of families 
and bifurcations) in comparison with the case of internal resonances. It is a fact that the 
complexity of the dynamics of continuation and evolution of families of periodic orbits 
depends on the existence of critical orbits of period T = 2kir, orbits of critical stability 
and asymmetric periodic orbits. For example the family S[ of the 2/1 internal resonance 
evolves smoothly as p changes and no any bifurcations are obtained in the domain (p^, 00). 
Also we can show that a similar smooth continuation and evolution holds for the family 
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7( e ) (parts J( ea ) and I^b) in Fig. As p starts from zero and tends to infinity, the family 
evolves smoothly and tends to the family (parts Jj a and Ij a ). No any bifurcations take 
place, while the collision orbit, which separates the parts (a) and (b) of the families, exists 
for all values of p. This family is equivalent to the family of (tt, tt)- corotations indicated 
by Beauge et al (2006) and has a similar evolution as that of the 3/1 family S2 of (tt, 0)- 
corotations shown in Voyatzis (2008). 

Apart from collision-bifurcations, structural changes in the characteristic curves of 
periodic orbits occur when their bifurcation points disappear. We showed how the asym- 
metric family A123, which starts and ends at bifurcation points on the symmetric family 
5*i, shrinks and finally disappears after the collapse of both bifurcation points. An in- 
teresting case occurs when the starting and the ending points of a family are bifurcation 
points that belong to different families. In this case it is possible, for a critical value of p, 
one of the points to disappear or to become a bifurcation point for another family. Such 
a bifurcation has been observed in the 3/1 resonance (Voyatzis, 2008) and cause global 
changes in the structure of the characteristic curves. 

Although we presented results that are associated with the 2/1 resonance, we can 
claim that similar features of continuation are present in other resonances. Our numerical 
study can not be extended efficiently up to high eccentricities (ei — > 1 or e 2 — > 1), due 
to computation restrictions. We believe that the features of the families in these regions 
can be revealed by considering the rectilinear restricted problem beside the circular and 
the elliptic one. 



References 

[1] Beauge, C, Ferraz-Mello, S. and Michtchenko, T.: 2003, 'Extrasolar Planets in Mean- 
Motion Resonance: Apses Alignment and Asymmetric Stationary Solutions', ApJ, 593, 
1124-1133. 

[2] Beauge, C, Ferraz-Mello, S., & Michtchenko, T.: 2006, Planetary migration and 
extrasolar planets in the 2/1 mean-motion resoance, MNRAS, 365, 1160-1170 

[3] Bozis, G. and Hadjidemetriou J.D. : 1976, On the continution of periodic orbits from 
the restricted to the general three body problem, Cel.Mech., 13, 127-136 

[4] Broucke R. A. : 1969a, Periodic Orbits in the Elliptic restricted three-body problem, 
JPL Tech.Rep., No 32-1360. 

[5] Broucke R. A. : 1969b, Stability of Periodic Orbits in the Elliptic restricted three-body 
problem, AIAA 7, 1003-1009. 

[6] Bruno, A.D. : 1994, The restricted Three body problem: Plane periodic orbits, Walter 
de Guyter, Berlin-New York. 



17 



[7] Bruno, A.D. and Varin, V.P.: 2006, On families of periodic solutions in the restricted 
TBP, Cel.Mech.Dyn.Astr., 95, 27-54. 

[8] Bruno, A.D. and Varin, V.P.: 2007, Periodic solutions of the restricted TBP for small 
mass ratio, Applied Math.&Mech., 71, 933-960. 

[9] Ferraz-Mello, S., Beauge, C. & Michtchenko, T.: 2003, Evolution of migrating planets 
in resonance, Cel.Mech.Dyn.Astr., 87, 99 

[10] Ferraz-Mello, S., Michtchenko, T. and Beauge, C: 2005, 'The orbits of the extrasolar 
planets HD 82943c and b', ApJ, 621, 473-481. 

[11] Hadjidemetriou, J.D. : 1975, 'The continuation of periodic orbits from the restricted 
to the general three-body problem', Cel.Mech.Dyn.Astr., 12, 155-174. 

[12] Hadjidemetriou J.D. and Christides, T. : 1975, Families of periodic orbits in the 
planar three body problem, Cel.Mech., 12, 175-187. 

[13] Hadjidemetriou, J.D. : 1976, 'Families of Periodic Planetary Type Orbits in the 
three-body Problem and their stability', Astrophys.Sp.Sci., 40, 201-224. 

[14] Hadjidemetriou, J.D.: 1993, Resonant motion in the restricted three body problem. 
Cel.Mech.Dyn.Astr. 56, 201-219 

[15] Hadjidemetriou, J.D. : 2002, 'Resonant periodic motion and the stability of extrasolar 
planetary systems', Cel.Mech.Dyn.Astr., 83, 141-154. 

[16] Hadjidemetriou, J.D., 2006: Symmetric and asymmetric librations in extrasolar plan- 
etary systems: a global view, Cel.Mech.Dyn.Astr., 95 , 225-244 

[17] Haghighipour, N., Couetdic, J., Varadi, F., Moore, W.B.: 2003. Stable 1:2 resonant 
periodic orbits in elliptic three-body systems. ApJ 596, 1332-1340. 

[18] Henon, M : 1997, Generating families in the restricted three body problem, LPN, 
m52, 1-278 

[19] Ji, J., Liu, L., Kinoshita, H., Zhou, J., Nakai, H. and Li, G. : 2003 The librat- 
ing companions in HD37124, HD12661, HD82943, 47UM and GJ876: alignment or 
antialignment, ApJ, 591, L57-L60. 

[20] Lissauer, J.J. and Rivera, E.J. : 2001, 'Stability analysis of the Planetary System 
orbiting Andromedae', ApJ, 554, 1141-1150. 

[21] Maciejewski, A.J. and Rybicki, S.M. : 2004, Global bifurcations of periodic solutions 
of the restricted TBP, Cel.Mech.Dyn.Astr., 88, 293-324. 



18 



[22] Michtchenko, T. A., Beauge, C, Ferraz-Mello, S.: 2006, Stationary solutions in 
resonant extrasolar planetary systems, Cel.Mech.Dyn.Astr., 94 , 411 

[23] Papadakis, K. and Goudas, C. : 2006, Restricted three body problem: An Approxi- 
mation of its general solution, Astrophys.Sp.Sci., 305, 99-124. 

[24] Psychoyos, D. and Hadjidemetriou, J.D. : 2005, 'Dynamics of 2/1 resonant extrasolar 
systems. Application to HD82943 and Gliese876', Cel.Mech.Dyn.Astr., 92, 135-156 

[25] Tsiganis, KL, Varvoglis, H., and Hadjidemetriou, J.D. (2002a): Stable chaos in 
Higher-order Jovian Resonances, Icarus 155, 454-474. 

[26] Tsiganis, KL, Varvoglis, H., and Hadjidemetriou, J.D. (2002b): Stable chaos versus 
Kirkwood Gaps in the Asteroid Belt: A comparative Study of mean Motion Resonances, 
Icarus 159, 284-299. 

[27] Voyatzis, G., Kotoulas, T., Hadjidemetriou, J.D.: 2005, Symmetric and Non 
symmetric Periodic Orbits in the exterior Mean Motion Resonances with Neptune. 
Cel.Mech.Dyn.Astr. 91, 191-202. 

[28] Voyatzis, G., & Kotoulas T.A.: 2005, Planar periodic orbits in exterior resonances 
with Neptune, Planetary and Space Sci., 53, 1189-1199. 

[29] Voyatzis, G. and Hadjidemetriou , H.D. : 2005, 'Symmetric and asymmetric librations 
in planetary and satellite systems at the 2/1 resonance', Cel.Mech.Dyn.Astr. 93, 265- 
296. 

[30] Voyatzis, G. : 2008, Chaos, order and periodic orbits in 3:1 resonant planetary 
dynamics, ApJ, 675, 802-816. 

[31] Zhou, L.Y, Lehto, H.J., Sun, Y.S. & Zheng, J.Q. : 2004, Apsidal corotation in mean 
motion resoance: the 55 Cancri system as an example, MNRAS, 350, 1495-1502 



19 




Figure 5: Families of periodic orbits of the general problem (solid curves) for p = 0.01 in 
cases (a)-(c) (external resonance) and p = 100 in case (d) (internal resonance). The 
families of the circular restricted problem lie along the axis e\ — and e2 = for 
the external and internal resonance, respectively. The families of the elliptic restricted 
problem are indicated by the dashed curves. 
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Figure 6: The "collision-bifurcation" of the families A 2 and A% at p = p 2 and the formation 
of the new families A 2 % and A% 2 . a) p — 0.2 < pi b) p = 0.275 ~ p\ and c) p = 0.3 > p\. 
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Figure 7: The "collision-bifurcation" of the families A\ and A^ at p = p 2 and the 
formation of the new families A 4 and A 123 . a) p = 0.3 < p 2 b) p = 0.37 ~ p 2 and c) 
p = 0.4 > p 2 . 
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Figure 8: a) The evolution of the families A 12 3 and Si as p passes the critical value 
P3 = 1.034 and takes large values, b) The stability index 6 2 along the family S\, which 
determines the interval of instability (&2 < —2) and the bifurcation points (at 62 = ~ 2) 
of the family A123. 
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